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SUMMARY 


A  hollow  reflecting  cone  or  pyranid  can  be  used  to  concentrate  infrared 
radiation  on  to  a  snail  photo-detector.  This  note  describes  calculations  which 
predict  that  a  five  fold  increase  in  flux  density  should  be  readily  obtainable. 
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SUMMARY 


The  RSRE  Scanner  Optics  Ray-tracing  Programs  have  been  used  to  calculate  the 
transmittance  of  a  truncated  pyramid  reflector  that  would  increase  the  apparent 
size  of  an  infrared  detector  if  placed  immediately  over  it.  The  dependence  of 
flux  concentration  with  pyramid  geometry  and  ray  angles  is  presented  graphically 
and  it  is  shown  that  concentration  factors  greater  than  five  should  be  easy  to 
achieve . 

2  BACKGROUND 

The  use  of  a  hollow  cone  or  pyramid  with  reflecting  walls  to  concentrate 
radiation  on  to  a  small  area  detector  is  well  known  (ref.l).  Figure  1.  shows 
a  parallel  bundle  of  30  rays  incident  at  an  angle  of  10  degrees  on  to  a  3:1 
truncated  pyramid.  Only  8  rays  are  reflected  back  out  of  the  pyramid  and  so 
the  nett  ray  concentration  factor  is  (22/30)  x  3=2. 2:1  (rays/unit  length). 

In  the  quest  for  thermoelectrically  cooled  infrared  detectors  it  has  been 
determined  that  a  2  or  3  fold  increase  in  flux  density  at  the  detector  could 
dramatically  improve  the  chances  of  successful  operation  (ref. 2).  An  array  of 
pyramid  reflectors  could  be  made  to  overlay  a  detector  array  whose  inter-element 
spacings  were  increased  to  match  the  desired  concentration  ratio. 

This  memo  reports  how  the  transmittance  of  a  truncated  hollow  pyramid  was 
determined  by  ray-tracing  methods  for  several  pyramid  geometries  and  ray  angles. 

A  simple  model  was  used  to  describe  the  ray  angle  distribution  at  the  focus  of 
a  lens  and  numerical  integrations  were  carried  out  to  estimate  the  flux 
concentration  factors  at  the  chosen  truncation  ratios,  and  pyramid  side  angles 
for  several  lens  F-numbers.  Graphs  are  presented  (Figs. .4, 5)  which  show  how 
the  concentration  factor  is  affected  by  all  these  parameters. 

3  METHOD 

3.1  THE  RAY -TRACING 

A  square  pyramid  of  1/2  angle  6  lying  along  the  z-axis  with  its  apex  at  the 
origin  is  bounded  by  the  four  planes 

z  1  x/tan(0)  *  0  and  z  1  y/tan(0)  “  0 . (1) 

If  a  ray  with  direction  cosines  l,m,n  passes  through  a  point  (a,b,c)  then 
another  point  on  the  ray,  a  distance,  d,  from  (a,b,c)  is  at  (x,y,z)  where 

x  ■  a  ♦  d*l  y  •  b  ♦  d*m  z  ■  c  ♦  d*n  .  (2) 

Equations  (1)  and  (2)  can  be  combined  to  give  the  distances  along  the  ray 
from  (a,b,c)  to  the  four  planes  that  form  the  pyramid; 

d  ■  -(c*tan(0)  1  b)/(h*tan(0)  ±  m)  .  (3) 

The  HP9825  computer  subprogram  "FYR"  listed  in  the  Appendix  was  used  in  the 
RSRE  Ray-trace  Programs  (Ref .3)  to  locate  the  point  of  intersection  on  the 
appropriate  face  of  the  pyramid.  A  control  flag  was  set  if  the  intersection 
point  was  closer  to  the  apex  than  that  appropriate  to  the  chosen  truncation 
ratio.  A  different  control  flag  was  set  if  the  ray  direction  had  a  component 
in  the  +z  direction,  indicating  that  the  ray  had  failed  to  pass  through  the 
pyramid  and  was  travelling  away  from  the  detector. 
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The  procedure  "PYR"  outlined  above  can  be  used  to  find  out  if  a  given  ray 
ia  transmitted  or  reflected  by  a  truncated  square  reflecting  pyramid.  A 
ray  incident  on  the  open  base  of  the  pyramid  at  a  given  angle  can  be 
scanned  in  (say)  the  y-direction  to  locate  the  reflection/transmission 
boundaries  for  rays  of  that  angle  and  x-position.  Once  these  y-values 
have  been  located  for  several  x-positions  a  numerical  integration 
procedure  will  determine  vhat  proportion  of  a  parallel  ray-bundle  filling 
the  pyramid  base  will  pass  through  to  the  detector.  This  fraction,  F,  was 
calculated  by  the  program  "RT2#6"  listed  in  the  appendix  for  ray  angles  from 
0  to  30  degrees,  for  pyramid  angles  from  10  to  25  degrees,  and  for  trunca¬ 
tion  ratios  from  2:1  to  3.5:1.  The  results  are  plotted  as  F  vs.  Ray  angle 
in  Fig.2(A-D). 


3.3  FORMULAE  FOR  FLUX  DENSITY  IMPROVEMENT  FACTOR 


Let  a  circular  lens  aperture  be  uniformly  illuminated.  After  refraction 
through  the  lens  the  angle  between  a  ray  and  the  lens  axis  is  roughly 
proportional  to  the  ray  height  at  the  lens  aperture.  The  number  of  rays 
at  a  given  angle  is  therefore  approximately  proportional  to  the  ray  angle. 
If  the  maximum  ray  angle  is  a0  then  the  following  equations  can  be  written. 


Lens  F-number  ■  V(2sin(a0))  .  (4) 

«o 

Total  ray  flux  on  lens  ■  fLlada  ■  K*ae  .  (5) 

«'o  a 

Ray  flux  passing  through  pyramid  »  K/aF(a)da .  (6) 

Jo 

Dividing  eq(6)  by  (5)  will  determine  the  pyramid  transmittance  for  a 
focussed  beam  of  given  numerical  aperture  (NA  *  sin(a0)). 


ie.  T  -  (2/0  /oF(a)da  .  (7) 

J  o 


If  all  radiation  incident  on  the  base  of  a  pyramid  with  (linear)  truncation 
ratio  R  passes  out  of  the  apex  end  then  the  flux  density  would  be  increased 
by  R*  Taking  equation  (7)  into  account  the  nett  flux  density  increase 
due  to  a  hollow  pyramid  concentrator  is  therefore  given  by 

«  “o 
2  2  f 

Flux  Density  Concentration  ■  (2*R  /a0  I aF(a)da  .  (8) 

o 

4  DISCUSSION  OF  RESULTS 

Equations  (7)  and  (8)  have  been  evaluated  by  numerical  integration  and  the 
results  are  shown  graphically  in  Fig. 3  and  Fig. 4.  In  these  figures  the  pyramid 
geometry  is  characterised  by  prism  angle  and  truncation  ratio.  An  alternative 
parameter  is  the  pyramid  height,  expressed  in  units  of  detector  size.  Different 
values  of  angle  and  truncation  ratio  can  combine  to  give  the  pyramid  height  (h) 
since  h*(R-l)/(2*tan(f)) .  Flux  density  concentrations  factors  have  been 
extracted  from  fig. 4  and  re-plotted  against  truncation  ratio  for  certain  F- 
numbers  and  pyramid  heights.  Figure  5  clearly  shows  that  concentration  factors 
in  excess  of  5  should  be  obtainable  with  reasonable  pyramid  geometries. 


It  is  implied  in  eq.(6)  that  F(a)  depends  only  on  the  ray  incidence  angle  and  not 
on  its  azimuth  (eg.  horizontal,  vertical,  or  diagonal).  A  square  pyramid  of 
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course  has  4-fold  synmetry  about  its  axis  and  so  the  integration  in  eq.(6)  is 
not  strictly  accurate.  However  the  variation  of  F(o)  with  azimuth  for  a  given 
ray  angle  is  expected  to  be  small,  making  eq.(6)  and  the  others  derived  from  it, 
good  approximations. 

A  ray-tracing  model  could  be  established  for  a  particular  lens,  pyramid,  and 
detector  arrangement  to  confirm  the  expected  flux  density  improvements.  Such 
a  computation  would  launch  a  large  bundle  of  rays  through  the  lens  representing 
uniform  illumination  (in  position  and  direction)  and  counting  those  that  pass 
into  and  out  of  the  pyramid  concentrator.  This  was  done  for  an  F-1.5  Germanium 
doublet  lens  feeding  a  25  micron  detector  via  a  2.5:1  pyramid  reflector  100 
microns  long.  The  incident  ray  bundle  comprised  nearly  50,000  rays  (the  ray 
bundle  was  intentionally  made  to  spread  well  beyond  the  base  of  the  pyramid). 

62Z  of  the  rays  passed  into  the  pyramid,  and  93. 1Z  of  these  reached  the  detector 
plane.  This  implies  a  flux  density  concentration  of  ,931*(2.5) *2»5.8,  within 
4Z  of  the  estimate  of  5.6  from  fig.5B. 
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0:  "RT2*6":dsp  "SPECIAL  PYRAMID  TRACE  ( 15/12/B3) " ;oto  "PROGRAM" istp 


61:  "PYR "  :  cfo  8,9;ina  Bifor  I*-l  to  3;for  J=1  to  3  ;R  Tpl  ,  I*3J+7]-»Z 
62:  Bm+Z<AIJ1-Rlpl  ;  BI 1+3  l*7At  J+3  l-*pr  I+3I  .next  J-.next  I 

63:  if  BI6J>0;sfo  B:ret 

64:  ton  <  R  fpl  ,221  >-*Z  :ZBt  3  ]  -*r3  ;  ZB  [  61-*r4  ;sf  q  14 
65:  -<r3+BIl]  )/<r4+Bt41  )-»r5;  if  prnd  (r5  .  -8)  <  =  0  :  9e99->r5 
66:  -<r3-Bt  l])/<r4-Bf4])-»r6;if  prnd  <r6  ,-8)  <  =  0  ;9e99-»r6 
67  :  -( r 3+BI21  )/<r  4+BI5  J  )-»r7 :  if  pr  nd(r  7 , -8 )  <  =  0 , 9»99-»r7 
68:  -<r3-B£21  )/<r4-Bf  53  >-»r8  ;  if  prnd <r8  ,-8)  <  =  0  ;9e99-»r8 
69:  cfg  1 4 t win < r 5  . r6 , r7 .rB > 4r2 
70:  if  r2=9e99; sfg  8; ret 

71:  for  1  =  1  to  3  :  Bill  +r2Bt  1+3  l-»BII]-*r  ( 1+6)  :nex  t  I 

72:  if  abs<B£3])<U/tan<R£l .22] > ,sfq  9;ret 

73:  prnd  <BI31  +BI 1 1/Z . -6)  -»r 5 

74:  prnd<B£3]-Bm/Z,-6)-»r6 

75:  prnd<BI31+B[2]/Z.-6)->r7 

76:  prnd<B£3]-Bf23/Z,-6)4r8 

77:  <r5=0)+2(r6*0)+3(r7«0  )+4<r8=0  )-»r0 

78:  <r0=2)-<r0=l  )-*r7:(r8=4>-(r0=3)-»r8 

79:  -Z-»r9.:l/f<r7r7+r8r8+r9r9>42  :Zr7-»r7:ZrB-»r8:Zr9->r9 

80:  ina  C  ;  0-*r  1 0-»r  1 1 -»r  12 

81:  for  1  =  1  to  3:for  J=1  to  3  iRtpl  .7+3I  +  JJ-»Z 
82:  Cm+7B[J]-»cm  ;  r  <  9+1  )+Zr  (6+J) -»r  (9+1 )  -.next  J 

83:  CmtR[pl,IHCm:ntxt  I;oto  “REFLECT  "  :  if  abo  <N)#9e99  ;ot  o  “REFRACT" 


156:  "PROGRAM": 

157:  if  not  flal2;sfa  12;di«  Xf 0 : 20 1 , Y f 0 : 20 ] .FT  0 : 20 ] 

158:  ina  R  :  90-*R  1 1  , 4  ]4R  [1  . 51-»RI1  . 8  ]  ;9e99-»R  [  l  .20  ]  :-9e99-»R  [  l  .  9] 

159:  ell  'ShG'(l) 

160  :  32-*UI  1  ]  •,  1 6-»U  £21  ;B-»U[3]  ;4-»Ut4]  ;2-»U[51  ;l-»ll[6]  ;  .5-*U[7]  ;  .  25-411  r  B 1 
161:  "NEXT": 74Q 
162:  for  U=1  to  5 

163:  <U=1 ) 4(3/3. 5)  <U=2) +  (3/2 .5)  <W=3)  +  (3/2)  <U=4>4<3/1  .5)  (U=5)-»U 

164:  fxd  4 ; pr t  “RATIO" . 3/U:spc 

165:  for  E*10  to  25  bv  5;E-*RI1,221 

166:  pr t  "P=",E:spc  :04O 

167:  for  A*0  to  30  bv  2;in#  F 

168:  fxd  6 ;pr  t  "A«" ,A 

169:  for  P=1  to  -1  bv  -2 

170:  -cos<A)4Dt5]  :P«iin(A)4DI6]  t  3/tan  (R[l  .22])  4DC2] 

171:  . 14H ; f or  G=0  to  10;.3G4Dfl] 

172:  04B4C:644F;for  D=1  to  8 :3F/1284DE3] 

173:  ara  D4C:for  V=1  to  lOOO  iara  C4A-,cll  'RAYTRACF'  <  1 ) 

174:  if  flg8:F-UtDJ4Fjato  43 
175:  if  f]q9;F4UtD]4F;qto  42 
176:  next  V 
177:  next  D 

178:  DI33/64FIG41  ]4FIG4i  ] -.next  G-.next  P 
179:  fxd  6;ll4N:cll  'SIMP'jprt  "S=",S:spc 
180:  A4X1041401  -.S4YIO]  -.next  A 

181:  pr t  " - " 

182:  E4XI0] ;3/U4YlO] ircf  0+14Q , X E » 1  , V t » 1 
183:  next  E;next  U 
184:  end 


0:  “RT7930  PYRAMID  RFSUI.TS  WITH  SPLINfi  FIT  < 06/01 /B4 o t o  "P" 

1  :  "DER INT  “  :  cto  3 ;  1-*I 

2:  if  < I*HI  )  >N-1  sot  o  *3 

3:  .5<XII3-Xtl-1  ]-»X)/<XtI  +  l 1-Xri-1  3-*H)-*Pf  I] 

4:  2<  <  < Yl 1*1 l-Yt 1 1 )/(Xt 1*1 1-XI I ] >-< Yl 1 1-Yt  T  —  1  3 )  /X )/H*T  >  *St 1 1  :3T-»GI  7  1  :ot  o 
5:  0*Sm*StNl  ;B-4F3*U 
6:  0*U:2-*I 

7:  ut-sm-umsu-i  i-<  .5-Brii>sri*i]*Gm>*T 
6:  if  <obs<T)*H)  >U: H*U 

9:  sm+T-»sm 

10  :  if  ION-1  ;I*1-»I  :pto  -3 
11:  if  U>*E ;qto  -5 
12:  0*1 

13:  if  <I*1*I>  >N-1  ;gto  *2 

14:  <S[I  +  13-Sm>/<X[I41  l-XII])-»GtI]  ;oto  -1 
15:  if  N=0 ; gtc  +13 

16:  04 J 

17:  if  <  J*1*J>  >M;qto  +11 

IB:  1*1  :  if  < T t J]*T ) >*X [ 11  :ot o  ♦? 

19:  *fa  3 ; dsp  “ARG  OUT  OF  BOUND" ;ret 
20:  if  <1+1*1 )  >N;ato  -1 
21 :  if  T>XII1 sgto  -1 
22:  1-1*1 

23:  Tt  Jl-Xtll+HjTt Jl-XtI+1  1*T;HT*X 

24:  Stll+HGIIMS 

25:  <1/6*Z> <StI]*StI+1  l+Sl+U 

26:  <  <YtI+l  l-Ym)/(X[I  +  l]-Xm)4M)H+Y[n*XU4RIJ] 

27:  W+<H+T>U*ZXG[I]*DfJ]  ;qto  -10 
28:  0+I+A 

29:  if  < 1*1+1 )  >N-1  •, ret 

30:  A+<  .5(XtI+n-X[I]-*H)<Y[I]+Yri  +  l  D-<l/24)Ht3(Srn+Sri  +  l  1)  >*A;ato  -1 
31  : 

32:  “P": 

33:  21+N;din  XtNl  ,  YfNl  ,SfNl  , Gf N-1  1  ,Uf?l  1  ,Vf  21  1 
34:  100-»h:di«  T(MI  .BtMJ  .DIM] 

35:  fnt  ;wrt  705 . "IP700 .000 , 71 00 , 1 0400 - 
36:  wrt  705, "VSIS" -.csit  1.3, 1.7. 1.5 
37:  7-»Q 

38:  “NEXT": for  Y*1  to  2 

39:  1  . 1+Z  tscl  0 , 30  .  -Z<2-Y  ) ,  ZY 

40:  fxd  0  \  x  a  x  0 ,5 ,0 ,30  ,l;f xd  l;yox  0  ,  .  1  , 0 , 1  , 2  vf  xd  6 

41 :  for  C*1  to  4 

42:  ldf  Q+i+o,ur*i  ,UI»1 

43:  if  C=1  ifxd  ljplt  25 , .95,1 ;  lbl  "TRUNC=" ,Vl 1  I ;f xd  6 
44:  for  1*1  to  20  ill!  1*1  l+XIIl  :  VI  1*11  +YII)  :ne«  t  I 

45:  c*iz  1 ,1 .6,3/2;wrt  705 , "SMX " ; f or  1  =  1  to  16-,plt  XT  II  ,Yf  1 1  . 1  :nex  t  I 
46:  wrt  705.  “SMM 

47:  for  1=0  to  99  ;30I/99*Tt  1  +  1 1 -.next  I 
48:  16+N;100+M;le-6+E;g*b  "DER INT " 

49:  for  1*1  to  100',plt  Tl  1 1 ,Bf 1 1 jnext  I;pen 


50:  next  C 

51  :  next  Y  --j 

52:  end  U 


"RT2931  PYRAMID  RESULTS  MITH  SPLINE  FIT  <ii>  <09/01/84) " :oto  “P “ 

"DER INT  "  i  cfg  3 : 1  -* I 
if  <I-H-»I>>N-1  ;qto  +3 

.5<xiii-xii-i  i-»x)/<xti+i  j-x[i-n-»H>-»itm 

2(((Y[IMJ-Ym)/(X(I+l]-Xm)-<YIIl-Y[I-l)>/X)/H-tT>-*Srn  ;3T-*Prn  ;  a  t  o 

0-»Stn-»SIN]  ;B-4f3-»U 

0-»U;2-»I 

ui-sm-imsi i-i  i-<  .s-sriDsii+i  i+gtiiwt 

if  <«bo<T)-»H)>U;H-*U 

sm+T+stn 

I:  if  I9N-1 sI+l-H  :oto  -3 
:  if  U>=E;qto  -5 
(:  0-*I 

I:  if  <  1*1  -*I >  >N— 1  ;qto  +2 

I:  <StI-Ml-Sm)/<X[I  +  U-Xm)-»G[I);oto  -1 
i:  if  M*0;gto  +13 
> :  0-»J 

if  <  J+1-*J)  >M;ato  +11 
I:  1-*I ;  if  <TIJ)-»T)>*Xm  ;oto  +2 
►:  »fq  3 ; dsp  "ARG  OUT  OF  BOUND" ;ret 
I:  if  <  I  +  HI ) >N ;gto  -1 
:  if  T>XIU  ;gto  -1 
f:  I-1-*I 

I:  TI JJ-XCU-»H;TI Jl-XfI+l  1-*T;HT-*X 
I:  sm+HGtl  )-*S 
i:  M/6-»ZXSm+S[I  +  l]+S)-*U 

» :  <  <yi  i+n-Ym)/<xii+u-xm  )-»w>H^Ym-»xu-»Bi jj 

y+<H+T)U+ZXG[IHOm  ;qto  -10 
l!  0-»I-»A 

if  <I  +  1-»I>  >N-1  ;  ret 

I:  A+<  .  5 (  X 1 1  +  1  ]-XtII**H)  (YIIJ+Ytl-tl  )  )-<l/24)Ht3<Sf  D+STI  +  l  1  )  >-*A;ato  -1 


:  ‘P  "  i 

!  21*»N;di«  XtN]  ,  YCNJ  .SINJ  ,GfN-l  ]  ,UI21 1  ,Vt21 1 
1 00-»Mjdim  TIM1 ,BIM ] .DIM] 
fnt  ;wr t  705, "IP700, 800 ,7100 ,10400" 
wrt  705, "VS15";csiz  1.3. 1.7, 1.5 
7-»0 

"NEXT": for  Y=1  to  2;15-»B 
l.l-»Z;*cl  0 ,30  , -ZB<2-Y)  ,ZBY 

fxd  0  ;  xox  0,5,0, 30, 1 jf xd  0  ;  y«x  0,1,0,15,5; fxd  6 

for  C*  1  to  4 

ldf  0+l-*Q,Ut«]  ,VI«) 

if  C=l;fxd  1  ; pit  25,B,l;lbl  "TRUNO * , Vtl] ;f xd  6 
for  1  =  1  to  20;uri+l  ]-»xm  ;V[I*1  )xm-*rm  snext  I 

csiz  1 ,1 .6,3/2;wrt  705 , "SMX" ;f or  1*1  to  16;plt  Xtl) ,Y[ IJ  ,  1 >next  I 
urt  705, "SM" 

for  1*0  to  99  ; 30 1 /99-»T [  I  + 1  ]  '.next  I 

16-tN  j  1  00-»M  j  le-6-*E ;  aob  "DER  I  NT  " 

for  1*1  to  lOO'.plt  Till  ,BT  I  ]  *,next  I;  pen 

next  C 

next  Y 

end 


0:  "RT2432  PYRAMID  RESULTS  WITH  SPLINE  FIT  (iij)  < 09/0 t /B4 > " : a t o  -P " 

1  ;  "DERINT"  :cfa  3;1+I 

2:  if  < I  +  l+I  >  >N-1  ;ato  +3 

3:  .5<xt n-xti-i  i+x>/<xti+n-xn-i  i+m+Brn 

4:  2<<<YII+ll-YIIJ)/<XII+ll-XlIJ>-<YtIl-Y[  1-11  )/X  )/H-*T  )  +SIIJ  :  3T+G  1 1 )  :oto 
5:  O+SIll+StN]  ;B-4f3+W 
6:  0+U-.2+I 

7:  U(-SII1-B[ I1SI 1-1 1-(  ,5-BtIl  >SI I+l l+Gfll )+T 
8:  if  <abs<T)-»H)  >U :H+U 

9:  sm+T+sm 

10  :  if  I#N-1  :I+1+I ; at o  -3 
11 :  if  U>*E;gto  -5 
12:  0-«I 

13:  if  (I+l+I)  >N-1  ,ato  +2 

14:  (SII+11-Sm>/(XIl4ll-Xm)-»Gm:oto  -1 
15:  if  M=  0 ; g t o  +13 
16:  0+J 

17:  if  ( J+l+J) >M  ;gto  +11 

18:  1+1 : if  <T[ Jl+T) >=X(11 joto  +2 

19:  sfq  3 ; dsp  "ARG  OUT  OF  BOUND" ;ret 

20:  if  <1+1+1 > >N;oto  -1 

21 :  if  T>XII1 ;qto  -1 

22:  I-l+I 

23:  T[J]-Xm+HiT[J]-XII+l  l+T-.HT+X 

24:  Sm+HGIIl+S 

25:  (1/6+ZXSm+SII  +  l  )+S)+U 

26:  <<YII  +  11-YtU)/<X[I  +  l]-XtTJ  )+W)H+YI  Il+XU+Bt  J 1 
27:  U+<H+T)U+ZXGt Il+Df  J] ;qto  -10 
28:  0+I+A 

29:  if  (I+l+I) >N-1 ;re t 

30  :  A+<  ,5<XtI+l  l-Xtll+HXYfll+Yl  1  +  11  )-<  1 /24  )Ht3  ( Bf  Il+S  r  1  +  1  ]  )  )+A;ato  -1 

31  : 

32:  "P": 

33:  100+N;din  XIN1  ,YTN1  ,StNl  ,GfN-l 1  ,Ul?l 1  ,V121 1 
34:  1 00+Midin  TIM1 ,BtMl .DIM! 

35:  fnt  iwrt  705, "IP700, 800, 7100, 10400" 

36:  wrt  705 , "VSl 5" icsi z  1.3, 1.7, 1.5 
37:  7+Q 

38:  "NEXT" :for  Y=1  to  2;1+B 
39:  l.l+Z:scl  0 ,30 ,-ZB<2-Y> ,ZBY 
40:  csiz  1 .3,1 .7,1 .5 

41:  fxd  0;xox  0 , 5 , 0 ,30 , 1 ; f xd  1;yax  0 , . 1 , 0 , 1 , 2 ;f xd  6 
42:  for  C=1  to  4 

43:  ldf  0+1+0, Ut*J,V[«J*,Um+R;Um+P 
44:  if  Oljfxd  ljplt  25,B,l:lbl  "TRUNC*"  ,  Vtl)  :f  xd  6 
45:  for  1*1  to  20 ;U 1 1+11+Xf II ; Vf I+l 3XII 1 +Y 1 1 3 ;nex t  I 
46:  dsp  "recalculating  100  XI«J.Y[«J" 

47:  for  1*0  to  99 ;30I/99+Tt I+l J ;next  I 
48:  1 00+M; 1 6+N; 1 e-6+E  jgsb  "DERINT" 

49:  rdn  XfMl ,YCMJ ;ara  T+X}«ra  B+Y;dsp 
50 :  rdn  UIN1 ,VIN] 

51:  0+Ul 1 1+VI 1 J ; 1+F ; for  D=2  to  30  bv  2;0+M;F+l+F 
52:  for  1*1  to  M  :  <  1-1  >D/<M-1  )  +TI I J  -.next  I 
53:  int <3 . 3D)  +  1 +N ; l#-6+E ; qsb  "DERINT" 

54:  dsp  D,2A/Dt2+Z : D+UtF  )  iZ+VtF J :nex  t  D 
55:  rdn  XI163 ,Y116I ;ara  U+X;ara  V+Y 

56:  csiz  I ,1 .6.3/2;wrt  705 , "SMX " jf or  1*1  to  16jplt  XIII , YI II , 1 :nex t  I 
r,7 .  ur  *  a  cm'1 

58:  for  1*0  to  99 ;30I/99+T[ I+l 1 ;next  I 
59:  16+N> 1 00+M ; le-6+E ;asb  "DERINT" 

60:  for  1*1  to  100;plt  T t I ] ,B f I ] jnex t  I;pen 
61 :  next  C 
62:  next  Y 
63:  end 


0:  "RT2633  PYRAMID  RESULTS  WITH  SPLINF.  FIT  <iv>  <  09/01  /B4 ot  o  "P" 

1  !  "DERINT":cfo  3:1+1 
2:  if  <1+1+1 > >N-1 :oto  *3 

3:  .5<Xm-XII-1  J+X)/<XII+1 3-XtI-l  HHWMIJ 

4:  2<<<YII+11-YIIJ>/<XII+11-X[I1>-<YIII-YI  I-1J  >/X>/H+T  )+Sf  II  :3T+Gf  I J  ;oto  -i 
5:  0+SllI+SlNl ;B-4F3+U 
6:  0+U:2+I 

7:  w<-sii]-nmsti-n-<  .5-bi  insn+n+cmi+T 

6:  if  <#bs<T)+H> >U:H+U 

9:  sm+T+sm 

10:  if  ION-1 :I+1+I  toto  -3 
11:  if  U  >=E :gto  -5 
12:  0+1 

13:  if  <I+1+I>  >N-1  ioto  +2 

14:  <StI+ll-SII))/<XtI+ll-X[IJ>+Gm;oto  -1 
15:  if  M=0  •, gto  +13 
16:  0+J 

17:  if  < J+l+J) >M;qto  +11 

18:  1+1 tif  <  T l Jl+T) >=X [ 1 J ;ot o  +2 

19:  sfq  3 ; dsp  "ARG  OUT  OF  BOUND" ;ret 

20:  if  <1+1+1 > >N;oto  -1 

21 :  if  T>XII1 ;qto  -1 

22:  I-l+I 

23:  Tl J 1-X I I J+H ;T l JI-Xl I+l 3+T-.HT  +  X 

24:  SIIl+HGm+S 

25:  <l/6+Z><Sm+SII  +  ll+S)+U 

26:  <  <YII  +  l]-Y[I))/<XtI+l 3-XII1 >+W>H+YlI3+XU+&[ J3 
27:  U+<H+T)U+ZXGm+Pr  JJ  ;qto  -10 
28:  0+I+A 

29:  if  (I+l+I) >N-1 iret 

30:  A+<  .5<XU  +  1 1-XII1+H)  lYtll+YII  +  l  3)-<l/24)Ht3(Sri3+Sr  J  +  l  ] )  >+A;qto  -1 
31  : 

32 :  "P " : 

33:  100+N;din  XINJ , YINJ ,SfNJ .G1N-1 J ,Uf 21 J.VI211 
34:  100+M:din  TIM1 . BtM ] . DIM ) 

35:  fnt  ;wrt  705 , " IP700 ,800 ,71 00 , 1 0400 " 

36:  wr t  705 , "VSl 5" : c*i z  1.3, 1.7. 1.5 
37:  7+0 

38:  "NEXT": for  Y=1  to  2;10+B 
39:  1 . 1 +Z : sc  1  1 ,5,-ZB(2-Y> ,ZBY 
40:  csiz  1 .3,1 .7,1.5 

41:  fxd  0;xax  0 , .5,1 ,5 ,2;f xd  0;y*x  1 , 1 , 0 , 1 0 , 2  ;f xd  6 

42:  for  C=1  to  4:rdn  U121 I ,U121J ,X120 3 .Y120 3 

43:  ldf  O+l+Q.Ut*!  ,Vl*3  jVtll+R  ,UtU+P 

44:  if  C~1 jf xd  Ijplt  4,9,l:lbl  "TRUNC= " . V113 : f xd  6 

45:  for  1=1  to  20  tUl  I  +  l  3+XI II  I+UXF  Il+Yl  13  -.next  I 

46:  dsp  "recalculatina  100  Xl*],Yt»J" 

47:  for  1=0  to  99  ;30I/99+Tl I  +  l 3 ;n#x t  I 
48:  1 00+M; 1 6+N : le-6+E ;gsb  "DERINT" 

49:  rdn  XIM3 ,Y[M3 ;<r«  T+X;oro  B+Y;dsp 
50:  rdn  UIN3 ,VIN3 

51:  0+UI1I+VI1 3 ;l+F;for  D=2  to  30  by  2;0+M;F+l+F 
52:  for  1=1  to  M; < 1-1 )D/<H-1 >+Tl I  3  t  nex  t  I 
53:  int  <3 . 3D3  +  1+N ;  le-6+E  ', qsb  "DERINT" 

54:  dsp  D ,2ARR/Dt2+Z iX[N1+UIF3jZ+VIF3 :next  D 

55:  rdn  XI 163 ,Y1163 ;f or  1  =  1  to  13  :  V117-I 3+Yl  1 3 ;  1  /?r,in  <UT  17-11  >+Xf  1 3  ;nrx  t  I 
56:  csiz  1 , 1 . 6,3/2swr t  705, "SMX" jf or  1=1  to  13>plt  XIII, YITI.lt next  I 
57:  wrt  705  "SM" 

58:  for  I-o’to  99;1+<X1133-1 3I/99+TII+1 3 ;next  I 
59:  1 3+N : 1 0 0+M ) 1 e-6+E ; qsb  "DERINT" 

60:  for  1=1  to  100;plt  Tl 1 3 ,Bf 1 3 jnext  Itpen 
61 :  next  C 
62:  next  Y 
63:  end 


0:  "RT2434  PYRAMID  RESULTS  UITH  SPLINE  FIT  <»)  < 1 2/0 1 /84> " : o t o  "P" 

1  :  "DERINT":cfa  3;1-»I 
2:  if  <1+ 1  -»I ) >N-1  ;ot o  +3 

3:  .5<X[I]-X[I-ll-»X)/(XII  +  ll-Xri-1  l-*H)-»F»ril 

4:  2<  <  <  Y 1 1*1  l-YtI])/(XII+n-Xm)-(Ym-VII-lI)/X)/H-»T)-»S[I]  :3T-*G[I  J  jot  a  - 
5:  0-»SI1  3  -»S  t  N  3  i B-4f 3-»W 
6:  0-MJ-.2-H 

7:  u<-sin-Bmsn-n-<  .5-B[ii>sti+n+Gm>-»T 
8:  if  <ebs(T)-*H>  >U  :H-»U 

9:  sm+T-»sm 

10:  if  I*N-1  ;I+l-*I;oto  -3 

11 :  if  U>=E ;gto  -5 

12:  0-»I 

13:  if  <I  +  1-»I>  >N-1  lato  *2 

14:  <StI  +  l]-Sm)/(XII  +  l]-XII))4Gm  ;oto  -1 
15:  if  M=0;gto  +13 
16:  0-*J 

17:  if  <J+l-»J)>M;gto  +11 

18:  1-»I  :  if  <TIJ1-»T)  >=Xt  1  J  lato  +2 

19:  sfa  3 ;  dsp  "ARG  OUT  OF  BOUND"  ; ret 

20:  if  (I  +  l-»I)>N;gto  -1 

21 :  if  T)X[I] jgto  -1 

22:  I-1-»I 

23:  TlJ]-Xm-»H;Tl  J3-XII+1 1-»T  ;HT-»X 

24:  Sm+HGIIHS 

25:  <i/6-*zxsm+sti  +  n+s>-»u 

26:  <  <YtI+l 3-YII3  >/<X[I+l 3-XJI3  YII3+XU-HH  J3 

27:  «+<H+T)U+ZXGm-tPtJ];flto  -10 
28:  0-»I-»A 

29:  if  <I  +  1-*I>  >N-1  ;ret 

30  :  A+<  .5<XII  +  1 3  -Xt  1 3 -»H  >  (YIIl+YII  +  1  ] )  -  ( 1  /24  >Ht3  (S  r  I3  +  SIJ  +  1  3  >>-»A;ato  -1 

31  : 

32:  "P": 

33:  1 00-»N ;  din  XfNl  ,  YTN3  .SIN3  ,GfN-l  3  ,UI21 3  ,Vt21 3 

34;  10  0-*M:din  TIM1  .BIM3  .DIM  3 

35:  din  R14.4.43 

36:  “NEXT " :  f  or  0=8  to  8+16 

37:  rdn  Ut21 3 ,V[21 3 .XI203 .YI203 

38:  ldf  O.UI«3  ,Vl*3  ;VI  1  3-»R:Utl  3-»P 

39:  for  1  =  1  to  20  ;  UT 1  +  1  34X1 1 3  ;  Vf  1+1  3XtI3-»Ym  -.next  I 
40:  dsp  "recalculating  100  Xt*3,YI*3" 

41:  for  1  =  0  to  99  ;30J/99-tTI  I  +  l  3  ;nex  t  I 
42:  100-»M;16-»N:  le-6-»E;osb  "DERINT" 

43:  rdn  XfHl  ,Y[M3  ;ara  T-»X;ara  B-*Y;dsp 
44:  rdn  UIN3 ,VIN3 

45:  0-»Utl  3-mi  3  ;  l-*F;for  D=2  to  30  bv  2 ;  0-*M  ;F+1 -»F 
46:  for  1  =  1  to  H;(I-l)D/(H-lMTII3:ne»t  I 
47:  int<3.3D)+l-*N;le-6-*E;qsb  "DERINT" 

48:  dsp  D  .2ARR/Dt2-*Z:XIN3-»UrF3  ;Z-*VtF3  -.next  D 

49:  rdn  X[  163  ,YU63  ;for  1=1  to  13 ;  VI 17-1 1-»YII3 ;  l/2sin  (Uri7-1 3  >-»X  III  ;nex  t  I 
50:  for  1=1  to  4 ;  !♦  .  5<  1-1  )-»TI  I  3  :nex  t  I 
51:  13-*N;4-»M;le-6-»E;osb  "DERINT" 

52:  for  1  =  1  to  4  :BI  I  i-*RI  1 ,2R-3, .  2P  -1  3  ;ne  x  t  I 
53:  next  O 
54:  end 
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"RT2436  PYRAMID  RESULTS  WITH  SPLINF  FIT  (mi)  ( 27/0 1 /B4 ) " ; o t o  "P " 
"DERINT"  :cfo  3-.1  +  I 
if  (I+l+I) >N-1 ioto  +3 
.5(xin-xn-i  i+x>/<xti+i3-xti-i  i+H>+pm 

2<  <  (YI  1*1  l-YIIl)/<XII  +  ll-XtIl)-(Ym-YCI-l  1>/X)/H+T)+SII )  ;3T+Gi  I )  oto 
O+Sll 1+SIN1  ; 8-4F3+W 
0+U;2+I 

W(-SIIl-BlI3SlI-ll-(  .5-Bl 1 1  >Sf  I  +  l  l+GTIll+T 
if  (abs(T)+H)>U:H+U 
sm+T+sm 

if  I6N-1  jl+l+l  -.oto  -3 
if  U>*F; oto  -5 
0+1 

if  (l+l+I) >N-1  into  +2 
(sii+ii-sm>/(xti+i  i-xiid+giii  -.oto  -1 
if  M=0  ;gto  +13 
0+J 

if  (J+l+J)>M;qto  +11 
1+1 :  if  (Tl Jl+T) >=Xt 1 J ;ot  o  +2 
sfa  3 ;  dsp  "ARG  OUT  OF  BOUND'*  jpet 
if  < I+l+I ) >N i oto  -1 
if  T>X(I1 ;qto  -1 
I-l+I 

T I J  3 -X  til +H ; T IJ1 -X 1 1  + 1 3  +T ; HT+X 
S 1 1 3+HG 1 1 3+S 
(1/6+2 )(StI3+StI+l 1+S3+U 

( (Yll+l ]-Y[I3)/(XtI+l 3-XlI3)+W)H+YII3+XU+BIJ3 
U+(H+T 3U+ZXGI I 3+DI J3 ;oto  -10 
0+I+A 

if  ( I+l+I ) >N-1 ;ret 

A+(  .5(X[I  +  1  3-XU3+H)  (YtI3+Yf  I  +  l  3  >- ( 1 /24 >Ht3 ( S Til +S r T+l  ]  )  )+A.ato  -1 

NRII  , 

1 00+N ; din  XIN3.YrNl.SfN3 .GIN-1  3 ,UI21 3 ,VI21 3 
1 00+M :din  TIM 3  .BIM  3 ,DIM1 
din  R  r  4 , 4 , 4  3 .014,41 
ldf  35.RI*3 

fnt  j ur t  705, "IP700, 800, 7100, 10400" 

Mr t  705, "VS5":csiz  1.3, 1.7, 1.5 

0+0 

"NEXT" : for  Y=1  to  2;10+B ,0+1+0 
1 . 1+Z; scl  2,3.5. -ZB( 2-Y> . ZBY 
csiz  i . 1 , 1 . 7, 1 . 5 

fxd  l;xax  0 , . 1  ,2 ,3 . 5 ,5 ; f xd  0;yax  2 , 1 , 0 , 1 0 , 2 ; f xd  6 
for  L=2  to  6  bv  2 
nax(2,2Ltan(10>+l >+rl+Tl 1 3 
nin(3.5,2Ltan(25)+l )+r2+Tt  4  3 
rl+(r2-rl 3/3+TI23 }rl+2(r2-rl >/3+TI33 
for  C=1  to  4 

for  1  =  1  to  4;RI0,I  .Cl+YIIl  ;  ,5(3+I)+xri3 ;npxt  I 
4+M+Nj 1 e-6+E ; osb  "DERINT" 
for  1=1  to  4 ;BI 1 3+DIC , 1 1 ;ncx t  Itnext  C 
for  C=1  to  4 

for  1  =  1  to  4  tOU  .C3+YII1  ;5(  I  +  l  )  +X 1 1 3  mrxt  I 
atn( . 5(TIC1-1 )/L)+P+Tll 1 
l+M;4+N;l*-6+E ;psb  "DERINT" 

BI 1 3+OI 1 , Cl ;next  Csrl+TI13 
for  1=1  to  4; 011  .I1+YII3  ; Tlll+Xfll  ;nnt  I 
csiz  .8,1. 6,3/2 
if  L=2 (Mr t  705, "SM2" 
if  L=4 ;Mr t  705, "8M4" 
if  L*6 imp t  705, "SM6" 

for  1  =  1  to  4 i p 1 1  XI 1 1  ,YI 1 3 , 1 ;n«x t  Iswrt  705. "SM" 
for  1=0  to  49:rl+(r2-r 1 1I/49+TII+1 3 :next  I 
50+M j 4+N ; l*-6+E ;qsb  "DERINT" 
for  1=1  to  50;plt  Til  1 ,BI 1 3 jnox t  I;pen 

66:  next  Ljnoxt  Y(*nd 
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